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ABSTRACT 

Aims. We investigate the influence of stellar bulges on the star formation and morphology of disc galaxies that suffer from ram 
pressure. Several tree-SPH (smoothed particle hydrodynamics) simulations have been carried out to study the dependence of the star 
formation rate on the mass and size of a stellar bulge. In addition, different strengths of ram pressure and different alignments of the 
disc with respect to the intra-cluster medium (ICM) are applied. 

Methods. The simulations were carried out with the combined N-body/hydrodynamic code GADGET-2 with radiative cooling and a 
recipe for star formation. The same galaxy with different bulge sizes was used to accomplish 31 simulations with varying inclination 
angles and surrounding gas densities of 10~ 27 gcirT 3 and 10~ 28 gcrrT 3 . For all the simulations a relative velocity of 1000 kms~' for 
the galaxies and an initial gas temperature for the ICM of 10 7 K were applied. Besides galaxies flying edge-on and face-on through 
the surrounding gas, various disc tilt angles in between were used. To allow a comparison, the galaxies with the different bulges were 
also evolved in isolation to contrast the star formation rates. Furthermore, the influence of different disc gas mass fractions has been 
investigated. 

Results. As claimed in previous works, when ram pressure is acting on a galaxy, the star formation rate (SFR) is enhanced and rises up 
to four times with increasing ICM density compared to galaxies that evolve in isolation. However, a bulge suppresses the SFR when 
the same ram pressure is applied. Consequently, fewer new stars are formed because the SFR can be lowered by up to 2MQvr _1 . 
Furthermore, the denser the surrounding gas, the more inter-stellar medium (ISM) is stripped. While at an ICM density of 10~ 28 g cirT 3 
about 30% of the ISM is stripped, the galaxy is almost completely (more than 90%) stripped when an ICM density of 10~ 27 gciir 3 is 
applied. But again, a bulge prevents the stripping of the ISM and reduces the amount being stripped by up to 10%. Thereby, fewer 
stars are formed in the wake if the galaxy contains a bulge. The dependence of the SFR on the disc tilt angle is not very pronounced. 
Hereby a slight trend of decreasing star formation with increasing inclination angle can be determined. Furthermore, with increasing 
disc tilt angles, less gas is stripped and therefore fewer stars are formed in the wake. Reducing the disc gas mass fraction results in 
a lower SFR when the galaxies evolve in vacuum. On the other hand, the enhancement of the SFR in case of acting ram pressure is 
less pronounced with increasing gas mass fraction. Moreover, the fractional amount of stripped gas does not depend on the gas mass 
fraction. 

Key words. Galaxies:clusters:general - Galaxies:ISM - Galaxies :bulges - Galaxies: spiral - Galaxies:clusters:intracluster medium - 
Methods:numerical 



1. Introduction 

As is known since the 1970s, the star formation rate of galaxies 
depends heavily on their environment. As stated by Butcher & 
Oemler ( 1978), cluster galaxies are found to be redder than field 
galaxies. This implies that cluster galaxies must form stars at a 
lower rate than their counterparts in the field. Although the star 
formation rate is lower in cluster galaxies, they also tend to have 
less gas than galaxies in less dense environments. 
Furthermore, an excess of blue objects has been found at 
redshift z - 0.5, compared to lower redshifts. Moreover, spiral 
galaxies have a higher star formation rate than ellipticals. 
Ellipticals on the other hand are more numerous in the centre of 
a galaxy cluster ( Barkhouse et al.|2009| l. Furthermore, Dressier 
|et al.| ( |1997| l found that between redshift z = 0.5 and z = 
the fraction of SO galaxies is strongly increasing. Hence, much 
more early-type galaxies can be found in clusters, which could 
be explained by different galaxy transformation processes in 
clusters of galaxies. 

Many different mechanisms have been studied in this context to 



explain the change in morphology and star formation of cluster 



galaxies. Among others, galaxy harassment (e.g. Moore et al. 
1998) and strangulation/starvation (e.g. Larson et al.||1980|l are 



discussed mechanisms. Furthermore, galaxy-galaxy interactions 
(e.g. |Sulentic|[T976l |Bushouse|[l987j |Solomon & liagel[l98"8"l 
Kewley et al.||2006 ) as well as ram-pressure stripping ( |Gunn| 
& Gott 1972) can be found. Because these different processes 
can happen at the same time, it is crucial to separate these 
processes which can be done by studying the evolution of the 
star formation rate in dependence on the environment (e.g. 
Poggianti et al.|2006) . 

Clusters of galaxies contain a hot (« 10 7 K) and very thin 
(10~ 29 - 10~ 26 gcm~ 3 ) gas, the so-called intra-cluster medium, 
detected by X-ray observations. Galaxies that are moving in the 
cluster's potential and consequently are exposed to the ICM, 
are feeling the ram pressure of this hot gas and the inter-stellar 
medium is stripped if the ram pressure exceeds the gravitational 
restoring force, according to the Gunn and Gott criterion ( |Gunn 
& Gott T972J. Additionally, ram-pressure stripping is also a 
mechanism for enriching the ICM with metals (e.g. Schindler 
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|& Diaferio|2008l |Lovisari et al.|2009| |Durret et al.|201 1) . The 
influence of ram-pressure stripping on the morphology and 
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star formation was exhaustively investigated in the last decade 



(e.g.|Gavazzi et al.|2001||Yoshida et al.|2004)|Scott et al.| 2010 



Kenney et al. 2004|i, by investigating HI deficiencies and Ha 



and CO observations to investigate regions of star formation. 
This showed that interactions with the ICM cause the galaxies to 
have morphological disturbances of the gas disc and extensive 
gas loss in the central part of a galaxy cluster (e.g. Boselli et al. 
2009 ). In addition, star formation has also been found outside 



the galaxies themselves (e.g. Hester et al. 2010). 
Another approach to study the different interaction processes 
of galaxies in a cluster are numerical simulations. Tidal in- 
teractions between galaxies and galaxy-galaxy mergers were 
investigated by Kronberger et al. ( |2006| [2007). Additionally, 
they provided studies of the influence on the rotation curves 



and velocity fields of spiral galaxies (Kronberger et al. 2007 
Kapf erer et al.| |2006| ). Furthermore, it has been found that 
nonaxisymmetric and nonbisymmetric features are introduced 
by tidal interactions (Ru bin et al.|1 999). Ram-pressure stripping 
simulations were also performed by several groups, which found 
a good agreement of the Gunn & Gott criterion in hydrodynamic 
simulations as well JAbadi et alT]|1999| IVollmer et al.1|2001[ 
Roediger & Hensler||2005| ). |Jachym et al .(fl2009| |2007|) investi- 



gated the influence of a time-varying ram pressure and different 
inclination angles of the galaxies relative to the movement 
direction of a galaxy. These simulations were performed mostly 
with a smoothed particle hydrodynamics (SPH) prescription. 
Other methods, including sticky p articles (e.g. |Vollmer et aT 



2001 



and Eulerian grid-codes (e.g. Roedi ger & Bruggen|2006 
2007| were also used to simulate ram-pressure stripping. 
Because one of the main aims of this paper is to study the 
evolution of the star formation rate depending on the environ- 
ment, the simulation code includes a star formation recipe. A 
high pressure exerted by the ICM to the ISM can enhance the 
star formation rate significantly, as e.g. found by Kronberger 
|et al.| ( |2008| ). Other simulations also confirm the enhancement 
of the star formation rate when ram pressure is applied. Bekki 



& Couch (2003) demonstrated that the pressure of the ICM 



can induce the collapse of molecular clouds and consequently 
trigger a burst of star formation. Also using SPH, their recipe for 
star formation employs the local dynamical timescale and the 
sound crossing time of the gas particles to transform them into 



stellar particles. Furthermore, Bekki & Couch (2011) used a 
Schmidt-law with exponent y - 1 .5 to control the rate at which 
new stars are produced. Vollmer et al. (2012) also use a model 
for star formation in their simulation prescription, adopting 
the sticky-particles approach for the hydrodynamics with ram 
pressure applied by an additional acceleration to the particles. 
Stars are being formed during cloud collisions, represented by 
tracer-particles with zero mass. These particles are then evolved 
passively. 

In our approach, star formation processes are not restricted to 
the galaxy or some parts of the galaxy itself but are separately 
calculated for each particle, considering temperature and den- 
sity in a self-consistent way throughout the whole simulation 
domain. In the same way, new stellar particles are spawned 
that conserve the phase-space properties of the gas particle, 
containing the mass of newly formed stars. Thus, enhancements 
in the star formation rate due to ram pressure are obtained 
self-consistently. Already adopted for ram-pressure stripping 
simulations in Kronberg er et al.| ([2008 ), we describe this recipe 
in more detail in section 12. II 

In this work, especially the influence of the combination of a 
galaxy's morphology and ram -pressure stripping on the star for- 
mation rate is investigated. In this context, the special case of the 



influence of a stellar bulge is studied. We use different strengths 
of ram pressure and different model galaxies to decode a bulge's 
influence on ram-pressure stripping. The paper is organised as 
follows: in Sect. 2 we present the numerical setup. In Sect. 3 the 
results are presented and are compared to observations in Sect. 4. 
Finally, in Sect. 5 a summary with the main conclusions is given. 



2. Numerical setup 

2.1. Simulation implementation 

All simulations in this paper were performed with the N- 
body/hydrodynamic code GADGET-2 developed by Volker 
Springel (Springel 2005 ). The N-body part treats the collision- 
less dynamics of the dark matter and the stellar component. The 
gravitational force computation is performed using direct sum- 
mation for short range forces and a sophisticated tree code (as 
first introduced by Barn es & Hut|19 86 ) and treePM code for the 
long range forces. The gas of the galaxies and the ICM is treated 
hydrodynamically via SPH (Gingold & Monaghan 1977 Lucy 
[T977| ). 

Furthermore, the version of GADGET-2 used here includes ra- 



diative cooling ( Katz et aL|[19 96), a recipe for star formation, 
stellar feedback and galactic winds as described in Springel & 
|HernquisT ( 2003 ). This so-called hybrid model for star formation 
was already applied in Kronberger et al. ( 2008 1 and Kapferer 



et al. (2009| |2010|). Because the star formation is of particular 



interest in this paper, we present the main properties of the star 
formation model in the following. 

The star formation in the ISM is spatially not resolved but is 
described in a statistical manner. The continuous matter distri- 
bution is represented by SPH particles with spatially averaged 
hydrodynamic quantities. The star formation is thus calculated 
for each such particle. It is assumed that stars form from cold gas 
clouds, which cannot be resolved in the simulations either. A cer- 
tain fraction of mass of a single SPH particle is then considered 
as cold gas clouds and the total gas density can be partitioned 
as p — p c + p/j, namely the cold gas clouds and the hot ambi- 
ent medium. From these cold gas clouds, the stars are formed 
on a characteristic timescale f*, whereas a certain mass fraction 
/? is again immediately released in the form of supernovae from 
massive stars with M > 8 M0. As in the other works mentioned 
before, we adopt p 1 — 0.1 according to a Salpeter IMF with a 
slope of - 1 .35 with lower and upper limit of 0. 1 Mq and 40 Mq 
respectively. Then, the simple model of star formation and hence 
the time evolution of the stellar density p* is 



dp* 
dt 



Pc oP c 



(1) 



In addition to mass release, energy is also injected into the 
ISM by supernovae. The ambient gas is heated by this feed- 
back energy and the cold gas clouds can be evaporated. On the 
other hand, gas cools through the implemented radiative cool- 
ing process. A minimum temperature of 10 4 K can be reached. 
Therefore, a redistribution of the two different gas phases takes 
place. The density distribution of hot and cold gas evolving in 
time can consequently be written as 



dpc Pc . Q Pc 1 f k < \ 

-77--- A P— + Kel (Ph,U h ) 

dt t+ t+ Uh — u c 
and 
dp;, 



dt f+ f* u h - u c 



A ne t (Ph< u h) 



(2) 



(3) 
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for hot (ph) and cold (p c ) gas density, respectively, with u being 
the internal energy, A the efficiency of the evaporation process 
(depending on the density) and A net the cooling function (Katz 
|et al.|[T996] >. / is used to distinguish between ordinary cooling 
(/ = 1) and thermal instability (/ = 0) when star formation 
is going on. The switch to thermal instability cooling is made 
using a simple density threshold p t h, motivated by observations 
(Kennicut 1989), namely / = 1 for densities higher than p± and 
f = otherwise. 

All assumptions lead to a self-regulating star formation because 
the growth of cold gas clouds is balanced by supernova feed- 
back. 



in model galaxy ' lnb' 
2.6M yr 1 . 



with a current star formation rate of about 



Finally, a prescription for galactic winds, according to Springel 
& Hernquist ( 2003 ) is included. It is not a physically motivated 
model, but is trying phenomenologically to reproduce winds ac- 
cording to observational findings according to Martin (1999). 
It is assumed that the disc loses a certain amount of gas M w 
due to winds, which is proportional to the star formation rate, 
M w = 77M* with 77 = 2, according to the observations, and M* 
being the star formation rate of long-lived stars, namely with 
mass M < 8 M©. Additionally, a fixed fraction of the supernova 
energy of x = 0.25 is included in the wind, which contributes to 
the escape velocity of the wind from the disc. 

2.2. initial conditions of the model galaxies 

The model galaxies were created with a disc galaxy generator 
developed by Volker Springel (details can be found in Springel 
et al. 2005) according to Mo et al. ( 1998). Three model galaxies 
with different bulges were created, all of them with a total 
mass of 1.08 x 10 12 M©. The gaseous disc has an initial mass 
of 2.38 x 10 9 M©. The bulge of the galaxies is modelled with 
collisionless particles such as the dark-matter halo and the 
stellar disc. Only the mass of the bulge is varied, the other mass 
fractions are unaltered to investigate the influence of the bulge. 
The halo circular velocity for all galaxies is 160kms _1 and the 
stellar disc-scale length varies from 2.8 to 3.3 kpc according 
to the different size of the stellar bulge. The mass resolution, 
number of particles and total mass used for the different model 
galaxies is shown in Table[T] Note that the additional mass of the 
bulge of galaxies 01b and 05b is removed from the dark-matter 
halo. Thus, all galaxies have the same total mass. Furthermore, 
galaxy lnb has no bulge. 

To avoid numerical artefacts, we evolved the model galaxies for 
2 Gyr in isolation before they were used for further simulations 
of ram-pressure stripping. Because the initial condition gener- 
ator simply generates a density and velocity distribution of the 
different components of a model galaxy, the evolution in vac- 
uum, i.e. no surrounding gas, is necessary to gain self-regulated 
star formation. Furthermore, the gas distribution has settled to 
an equilibrium distribution with a spiral structure after 1 Gyr 
and a ring-like structure after 2 Gyr of evolution, as shown 
in Fig. 1(a) for model galaxy lnb. The radial density profile 



1(b) again shows model 
ution. In blue and red we 
the hot (> 10 4 K) 



2.1 



and density distribution plot in Fig. 
galaxy lnb after 2 Gyr of isolated evo^ 
display the two gas phases (see Sect 

and cold (< 10 4 K) gas of the model galaxies. The several local 
maxima originate from the pattern of the spiral structure of the 
disc. Note that we consider only the gas particles instead of the 
smooth density distribution for the calculation of the density 
profiles. Combined with the considered annuli and the limited 
mass resolution of the particles, the density distribution seems 
to have empty areas at a radius of 15, 20, and 25 kpc. 
During 2 Gyr of evolution, 4.2 x 10 9 M© new stars have formed 





Galaxy ' lnb' 


Galaxy '01b' 


Galaxy '05b' 


# of particles 








DM halo 


3 x 10 5 


3 x 10 5 


3 x 10 5 


gaseous disc 


2x 10 s 


2 x 10 5 


2 x 10 s 


stellar disc 


2 x 10 5 


2 x 10 5 


2 x 10 5 


stellar bulge 





2x 10 4 


1 x 10 5 


mass resolution 






[Mq /particle 








DM halo 


3.615 x 10" 


3.612 x 10" 


3.602 x 10" 


gaseous disc 


1.191 x 10 4 


1.191 x 10 4 


1.191 x 10 4 


stellar disc 


3.572 x 10 4 


3.572 x 10 4 


3.572 x 10 4 


stellar bulge 




4.762 x 10 4 


3.810 x 10 4 


total mass 








[Mq] 








DM halo 


1.084 x 10 12 


1.083 x 10 12 


1.080 x 10 12 


gaseous disc 


2.381 x 10 9 


2.381 x 10 9 


2.381 x 10 9 


stellar disc 


7.143 x 10 9 


7.143 x 10 9 


7.143 x 10 9 


stellar bulge 





9.524 x 10 s 


3.809 x 10 9 



Table 1. Initial conditions of the three model galaxies with num- 
ber of particles, mass resolution, and total mass of the different 
galaxy components. The total mass of all three model galaxies is 
the same. When including the bulges in model galaxies 01b and 
05b, less mass is consigned to the dark-matter halo when initial 
conditions are generated. 



2.3. Simulating the ram pressure 

To simulate a galaxy cluster environment, the model galaxies 
were exposed after 2 Gyr of evolution to an external wind, 
simulating the intra-cluster medium (ICM). For that purpose, a 
cuboid of gas particles was established, as depicted in Fig. [2] In 
this case, contrary to other ram-pressure stripping simulations, 
the galaxies themselves are moving instead of the wind and the 
ICM is stationary, which leads to the same result. 
The simulation domain was filled with gas particles of constant 
mass to achieve the desired density of the ICM. To avoid numer- 
ical artefacts due to the SPH scheme, the mass resolution of the 
particles representing the ICM was the same as that of the parti- 
cles representing the ISM in the model galaxy. Because for this 
setup, GADGET-2 obliges cubic simulation domains, a box with 
850 kpc side length was used. Filling this cube with particles of 
the given mass resolution would result in a huge amount of gas 
particles for the ICM and an unreasonable computing time. To 
avoid this problem, only in the centre a cuboid with a size of 
150 x 150 x 850 kpc was filled with gas particles, which had the 
same resolution as the model galaxies. The surrounding of this 
cuboid consisted of a less resolved medium with a 100 times 
higher mass per particle, which kept the inner simulation do- 
main, through which the model galaxies fly, stable. Additionally, 
periodic boundaries for the whole simulation domain were used. 
Galaxies moving with a velocity of 1000 km s~' cross the whole 
simulation domain in about 0.83 Gyr. The ICM does not change 
its density during this time by more than 4%. 
In the simulations, the model galaxies are moving with a con- 
stant velocity in a constant ambient medium. Since the ICM den- 
sity is constant during the whole simulation, this can be thought 
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hot gas 

— old stellar disc 
- newly formed stars 




10 15 20 

r [kpc] 
(a) 




Fig. 1. (a) Radial density profile of model galaxy lnb after 2 
Gyr of evolution in the vacuum. We show two gas phases used 
in the simulation, named hot (ionized) and cold gas. In black, 
the old stellar disc is plotted and we show in yellow the newly 
formed stars within the 2 Gyr of evolution. Owing to the ring- 
like structure and the method of calculating the profile, it seems 
that some gas is missing at a radius of 15, 20, and 25 kpc. (b) 
Face-on view of the galaxy's gas component. Already after 2 
Gyr, overdensities in concentric rings are formed. Furthermore, 
a spiral structure can be seen. 



of as galaxies orbiting a galaxy cluster, perceiving the same ram 
pressure throughout. The different configurations of ICM densi- 
ties, model galaxies, and inclination angles are depicted in Table 
[2] For all simulations, a relative velocity of 1000 km s _1 was 



used. Since Kapferer et al. ( 2009 ) found that the star formation 
rate depends much more strongly on the density of the ICM than 
on the relative velocity of a galaxy, we used only one velocity 
for the ram-pressure stripping simulations in this investigation. 
Furthermore, the ICM has a temperature of 10 7 K. 



Fig. 2. "Wind tunnel" used to simulate the ram-pressure strip- 
ping. The low-resolution gas cube has a side length of 850 kpc. 
The high-resolution wind tunnel in the centre has a dimension of 
150 x 150x850 kpc. 



3. Results 

3.1. The influence of a bulge and the inclination angle on the 
gaseous disc 

First, we investigated the effect of the ram pressure on the 
gaseous disc. As expected and shown in recent simulations on 
ram-pressure stripping, at the beginning of the interaction of 
the galaxy with the ICM, the spiral arms are being compressed 
and the ISM density increases in the first 50 Myr of acting ram 
pressure. Moreover, the stripping of the gaseous disc begins in 
the outskirts of the galaxy. 

In Fig. [3] we show the amount of stripped gas for the different 
simulations, and plot the fraction of gas in the wake of the 
galaxy to the total ISM gas mass versus simulation time. 
Obviously, when galaxies are flying face-on through the ICM, 
the gaseous disc is almost completely stripped after 100 Myr in 
the simulations with a high-density ICM of 10~ 27 gcirr 3 , and 
only a small amount of dense gas remains in the disc in the 
very centre. Owing to the rotation of the galaxy, filaments are 
forming in the stripped gas wake from the former spiral arms of 
the galaxy. After 200 Myr, very dense gaseous knots are forming 
in the wake. This scenario is also depicted for model galaxy lnb 
in Fig. [4] where the gas surface-density is plotted. The gas tail 
containing these knots is 40 kpc long after 200 Myr. After 500 
Myr of evolution the gas knots are ranging even to a distance 
of more than 300 kpc behind the disc. These self-gravitating 
gas knots are stable structures throughout nearly the whole 
simulation time although the mass of the knots is decreasing 
by 50-70 % as depicted for one distinct gas knot in Fig. [8] 
Moreover, the effect of a slipstream caused by the remaining gas 
disc can be observed. The filamentary structure is still present, 
albeit weakly. The gas density in these knots is comparable to 
the gas in the central region of the galaxy, reaching a value of 
2.66 x 10" 24 gem- 3 . 

In the-low density ICM (10~ 28 g cirT 3 ) simulations, filaments in 
the gas tail can only be seen after 200 Myr, a small amount of 
gas is stripped (20-30 %) and the main gaseous disc remains 
intact. After 500 Myr of acting ram pressure the filamentary 
structure reshapes to a few gaseous knots formed in the wake. 
Yet, the most dense regions are located at the centre of the 
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galaxy lnb, ICM density: 10 27 gem 3 galaxy lnb, ICM density: 10 2 *gcm 3 




Fig. 3. Fraction of the ISM, present in the wake of the galaxies, compared to the total amount of ISM is plotted versus simulation 
time, (a) and (c): high ram pressure, (b) and (d): low ram pressure. In (a) and (b), the behaviour of the different model galaxies 
is depicted to show the influence of a stellar bulge. In (c) and (d) we show the model galaxy lnb, flying with different inclination 
angles with respect to the wind. 



galaxy. 

The different model galaxies with distinct bulge masses behave 
differently when exposed to ram-pressure stripping, as can 



be seen in Fig. 3(a) for a high-density and Fig. 3(b) for a 



low-density ICM. Increasing the mass of the bulge, less gas is 
stripped, though the difference in the mass fraction is only a 
few percent for high ram pressure acting on a galaxy. If a low 
ram pressure is acting on the galaxy, the difference is more 
pronounced. In this special case, again the stripping is almost 
completed after 200 Myr. But with the low ram pressure acting, 
some part of the gas is falling back onto the disc, hence there is 
a local minimum at 300 Myr of simulation time. However, more 
gas is being stripped afterwards, which amounts to a total of a 
third of the gas. 

The impact of the inclination angle of a galaxy on the amount 
of stripped gas also depends on the ICM density, as depicted 
in Fig. |3(c)| and |3(d)| for model galaxy lnb. As expected, 
with increasing inclination, a galaxy loses less gas. While for 
inclination angles below 45° the stripping process is merely 
delayed, in the extreme case of galaxies flying edge-on through 



the ICM, only half of the gas is stripped in comparison to the 
face-on case. This influence of the inclination angle is the same 
for all model galaxies. 



3.2. The influence of a bulge and the inclination angle on the 
star formation rate 

In Fig. [5] the star formation rate of the whole simulation domain, 
including the wake, is plotted versus the simulation time for dif- 
ferent scenarios. For comparison, the SFR for the three model 
galaxies in vacuum is plotte d in F ig. |5(a)"| 



If a high ram pressure (Fig. |5(b)| ICM density of 10 27 gem 3 



with a face-on orientation) is acting on the galaxies, the SFR 
is enhanced up to a factor of 4 compared to the same galax- 
ies evolving in the field. Furthermore, a bulge decreases the star 
formation rate. For model galaxy 05b, the mean SFR is even 3 
M© yr _1 lower than that of model galaxy lnb. 
When a low ram pressure is acting (Fig. 5(c) ICM density of 



10 28 g cm 3 , with a face-on orientation), the star formation rate 
can be doubled. However, only a slight dependence on the bulge 
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Star formation rate (vacuum) 



galaxy lnb 

— galaxy 01b 

— galaxy 05b 




j.O 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
time [Gyr] 

(a) 



SFR,, ICM density 10 g cm , galaxy orientation: face-on 



SFR, ICM density 10 & gcm 3 , galaxy orientation: face-on 





SFR, galaxy lnb, ICM density V)-^gcmT 



SFR, galaxy lnb, ICM density 10~ 28 g cm- 





Fig. 5. Evolving star formation rates during the simulations. In Fig. (a) the star formation rates of the three model galaxies evolving 
in vacuum are shown, which are at about l.SMoyr" 1 for all galaxies. In Figs, (b) and (c) the star formation rates of the galaxies 
undergoing ram-pressure stripping are plotted, with an ICM density of (b) l(T 28 gcirr 3 and (c) 10~ 27 gcirT 3 . Already with a low 
ram pressure, the star formation is enhanced by a factor of two while with a high ram pressure the enhancement can reach a factor of 
4 after the first interaction of the galaxies with the ICM. Clearly, a bulge suppresses the star formation rate. In Figs, (d) and (e) the 
star formation rate of model galaxy ' lnb' undergoing ram-pressure stripping, flying with different inclination angles with respect to 
the ICM is shown. The star formation is suppressed if the galaxy is flying edge-on in case of high ram pressure. When a low ram 
pressure is acting on the galaxies, the inclination angle has almost no influence, the star formation rate is not affected. 
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# 


galaxy 




ICM density [g cm 3 ] 


inclination angle [°] 


1 


Inb 









2 


Inb 




10 -28 





3 


Inb 




1Q -28 


22.5 


4 


Inb 




10 -28 


45 


5 


Inb 




10 -28 


67.5 


6 


Inb 




1Q -28 


90 


7 


Inb 




10~ 27 





8 


Inb 




IO" 27 


45 


9 


Inb 




10~ 27 


90 


10 


01b 




o 




11 


01b 




10 -28 





12 


01b 




1Q -28 


22.5 


13 


01b 




1Q -28 


45 


14 


01b 




JQ-28 


67.5 


15 


01b 




1Q-28 


90 


16 


01b 




IO" 27 





17 


01b 




IO" 27 


45 


18 


01b 




io~ 27 


90 


19 


05b 




o 




20 


05b 




10 -28 





21 


05b 




1Q -28 


22.5 


22 


05b 




10 -28 


45 


23 


05b 




1Q-28 


67.5 


24 


05b 




1Q -28 


90 


25 


05b 




IO" 27 





26 


05b 




IO" 27 


45 


27 


05b 




io- 27 


90 


28 


Inb (10.0 


%) 







29 


Inb (17.5 


%) 
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Table 2. Properties 
#28-31, the disc gas 
from the initial 25% 



of the different simulations. In simulations 
fraction of the model galaxies was changed 
to 17.5% and 10%, respectively. 



mass can be found, again showing a lower SFR with increasing 
bulge mass. But a periodic, sinusoidal change in the star forma- 
tion rate can be observed because of the properties of the star 
formation recipe: in a periodic manner, the gas cools down and 
is compressed, hence new stars are formed. Consequently, the 
surrounding of the newly formed stars is heated up again due to 
the stellar feedback, which also decreases the density and hence 
lowers the SFR. Then, this process starts again. 
When the inclination angle of a galaxy relative to the movement 
is changed, the influence on the SFR is not very pronounced. In 
Figs. 



5(d) 



and 5(e) the SFR of model galaxy Inb is shown for 
a high and low ram pressure, respectively, and with different in- 
clination angles. With a low-density ICM and hence a low ram 
pressure, the sinusoidal shape of the star formation rate can be 
seen again and only a slight trend of decreasing star formation 
with increasing inclination angle can be observed. Obviously, 
the highest SFR can be seen in the face-on case, which is almost 
0.5 M0 higher than when the galaxy is tilted. If the inclination 
angle is increased, this effect becomes less pronounced. If there 
is high ram pressure, as already seen in the face-on case, the 
SFR is more enhanced for all inclination angles than in the case 
of low ram pressure. Again, the highest SFR can be seen in the 
face-on case, with decreasing SFR as the inclination angle in- 
creases. Interestingly, for a 45° angle compared to the edge-on 
case, the SFR is higher in the beginning of the simulation, while 
it is lower at the end, resulting in almost the same amount of new 
stars being produced. 



40 kpc 



Fig. 4. Surface density plot of the gaseous component of model 
galaxy Inb after 200 Myr of acting ram pressure. The filamen- 
tary and knot-like structure in the wake, which already has a 
length of about 40 kpc, is obvious. The galaxy is located on the 
left-hand side where only a bunch of dense gas remains in the 
central region. In the picture, the galaxy moves towards the left. 



3.3. The influence of the disc gas fraction on star formation 
rate and ram-pressure stripping 

In addition to considering the differences in SFR of a galaxy 
and the amount of stripped gas in a ram-pressure stripping 
scenario, the influence of the gas mass fraction on the disc needs 
to be scrutinised as well, which we did for model galaxy Inb. 
In Fig. |6| the evolving SFRs of simulations #1 and #28-30 are 
displayed. In Fig. |6(a)| the SFRs of model galaxy Inb evolving 
in vacuum are shown, with gas mass fractions of 25%, 17.5%, 



and 10% respectively. In Fig. 6(b) the SFRs of the same model 



galaxies undergoing ram-pressure stripping in an ICM with 
density 10~ 2c gcrrT 3 are displayed. Obviously, increasing the 
gas mass fraction in the disc leads to a higher star formation rate 
both when galaxies evolve in vacuum and when they undergo 
ram-pressure stripping. But with decreasing gas mass fraction, 
the SFR is more enhanced by ram-pressure. 
The fractional amount of stripped gas is independent of the 
gas mass fraction in the disc. Because the whole disc mass, 
including the stars, remains the same for the model galaxies, 
the restoring force according to the Gunn & Gott criterion also 
stays the same. Hence, these statements are also valid for the 
model galaxies with different bulge sizes. 



3.4. Distribution of the ISM and newly formed stars 

As shown before, when the ram pressure acting on a galaxy 
is high, almost all of its gas can be stripped. Nevertheless, the 
system of galaxy plus stripped gas still shows a high SFR. 
Obviously, the stars are formed in the stripped gas wake of the 
galaxy. 

In Fig. |7j four different timesteps of simulation #7 are shown 
with the gas and newly formed stars' surface density. Evidently, 
after 200 Myr, several dense gas knots are present in the wake 
with aggregated stars that were formed in these dense gas knots. 
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SFR, galaxy lnb, vacuum, different gas fractions 




(a) 



SFR, galaxy lnb, ICM density 10 28 p cm 3 , different gas fractions 

5r 




0.3 0.4 0.5 
time [Gyr] 

(b) 



Fig. 6. Star formation rates in simulations #1 and #28-30 are 
plotted. In (a), SFR of galaxy lnb with different gas fractions in 
the disc are plotted, evolving in vacuum. In (b), SFRs of the same 
three model galaxies undergoing ram-pressure stripping are dis- 
played. As expected, a lower gas mass fraction in the disc im- 
plies a lower star formation rate. In the ram-pressure stripping 
case, the star formation is enhanced in all cases, but with in- 
creasing gas mass fraction, the effect is less pronounced. 



Since the stars do not feel the ram pressure any more and are 
attracted by the galaxy's potential, they are falling towards and 
through the disc. Hence, stars are present also in front of the disc, 
as depicted in the timestep at 200 Myr. 

The first gas knots are formed after about 150 Myr of interaction 
with a dense ICM or rather when a high ram pressure is acting 
on a galaxy. The knots are present throughout the whole simu- 
lation for a total interaction time of the galaxy with the ICM of 
800 Myr. Nevertheless, the gas knots are losing mass, as shown 
in Fig. 8(a) for a gas knot taken from simulation #7. The mass 
loss is only due to star formation, thus the lost gas mass is trans- 
formed into new stars. The stellar feedback is not able to destroy 
the gas knots or to contribute to a significant outflow of the gas 
from the knots. Hence, cooling processes, self gravity and the 
external pressure from the ICM keep the gas knots stable over 
time and lead to an SFR within these knots throughout the whole 
simulation in the range of 0.025 to 0.01 M0yr _1 . In Fig. 8(b) 



the mass of newly formed stars in the distinct gas knot is shown 
throughout the simulation. Up to 10 6 Mq of newly formed stars 
can exist within the gas knot. Continuously, stars are falling back 
onto the galactic disc, creating streams of stars in the wake of the 
galaxy, as shown in timestep t = 200 Myr in Fig. [7] 
The amount of stars formed in the wake is a direct consequence 
of the amount of stripped gas, which differs for distinct bulges 
and different inclination angles of a galaxy and of course for dif- 
ferent ram pressure. With a low ram pressure, very few stars with 
a cumulative mass of « 4 x 10 M0 are present in the wake after 
500 Myr of interaction with the ICM. With increasing inclina- 
tion angle of a galaxy, also fewer stars are present in the wake as 
less gas is stripped. Consequently, also the presence of a bulge 
leads to less star formation in the wake, as expected. 
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Fig. 8. Mass of a distinctly chosen gas knot in the wake of ram- 
pressure stripped galaxy lnb in simulation #7 is shown with re- 
spect to simulation time. The gas clump is formed after 160 Myr 
of the onset of the ram pressure. Until the end of the simulation, 
the gas knot loses half of its mass due to star formation as shown 
in (a). The mass of newly formed stars, present in the gas knot, 
is displayed in (b). 
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4. Comparison with observations and other 
ram-pressure stripping simulations 



In this work, an SPH scheme as in Kronberger et al. (2006 



2007) and Kapferer et al. (2009) was used. In comparison to 
other ram-pressure stripping simulations that use a grid code 
(Roediger & Briiggen 2006, 2007), the behaviour is quite dif- 
ferent. In the grid codes, Kelvin-Helmholtz or Rayleigh-Taylor 
instabilities are much better resolved. In this case, fewer or even 
no gas clumps are formed. Different approaches for solving the 
hydrodynamic equations such as modified SPH schemes (e.g. 
Hefi & Springel|[20T0) |Read et al.||20l0|) or hybrid eulerian- 



lagrangian codes as in Springel (2010) also show less dumpiness 
in a galaxy's wake. However, self-gravity, cooling processes, and 
the external pressure from the ICM can indeed contribute to the 
stability of the gas clumps in a galaxy's wake, which then lead 
to star formation occurring outside of galaxies. Furthermore, the 
initial spiral structure of the gaseous component in the model 
galaxies, which are in a first moment compressed, lead to the fil- 
amentary structure and the final formation of the gas clumps. 
Galaxies undergoing ram-pressure stripping, as shown in our 
simulations, have been observed e.g. in the Virgo Cluster. The 
VIVA survey (VLA imaging of atomic gas in the Virgo Cluster) 



shows galaxies NGC440 2 and NGC4522 ( |Crowl et al 
Crowl & KenrTey] |2006]l 



2005 



es for 



which are prominent exampl 
interactions of the interstellar with the intra-cluster medium. 
For these cases, Crowl and collaborators showed the possibil- 
ity that galaxies can lose almost all of their gas in the case of 
high ram-pressure, as showed also apparently in the simulations. 
Furthermore, it can be seen that ram-pressure stripping can trig- 
ger a starburst with a subsequent cessation of the SFR, especially 
if the galaxy loses all of its gas. Such a strong ram-pressure ef- 
fect can not only be found in central regions of a galaxy cluster 
but also in regions well outside the cluster core, which affirms 
the validity of the orbiting scenario in our simulations. 
Another aspect is the transformation of spiral galaxies into SO or 
early-type galaxies. Tra n et ah] ( [20 05 ) found by studying galax- 
ies in MS2053 at redshift z = 0.6 that the Butcher-Oemler effect 
is linked to galaxy infall. In this case, as expected and seen in 
the simulations, ram-pressure stripping also plays an important 
role. Furthermore, [Wolfe t al. (2009) found in the context of the 
STAGES survey that the star formation rate in the observed clus- 
ter is clearly lower than in the field. Blue spirals in the field might 
be transformed first into red spirals and eventually into SOs in the 
cluster centre during the infall. Again, ram-pressure stripping, as 
shown in the simulations, might play an important role. 
The influence of the inclination angle of a galaxy's movement 
with respect to the ICM on the amount of stripped gas, which 
we have found, is confirmed by B osch et aT] ( |2012[ ). The galaxies 
in the Abell 901/902 supercluster show a clear correlation in the 
sense that a more face-on infall of a galaxy leads to a smaller gas 
disc with respect to the stellar disc, indicating a higher amount 
of stripped gas. 

Another prominent example is the galaxy IC 3418 in the Virgo 
cluster, observed with GALEX (Marti net al.|2005| ). It has been 
found that this galaxy has an ultraviolet tail of bright knots and 
diffuse emission with a length of about 17 kpc. Near ultraviolet, 
far ultraviolet and Ha emissions indicate that star formation is 
ongoing ( Gavazzi et al.||2006 i. Hester et al. ( 2010[ ) suggested, 
that IC 3418 must have been ram-pressure stripped in order to 
lose that amount of gas. As in the simulations, these authors 
claim that the gas clumps must have been formed in the stripped 
wake, aided by turbulence. Furthermore, the suggested star for- 
mation is taking place in these dense gas knots as well. 



More evidence for gas knots forming in the wake of ram- 
pressure stripped galaxies is provided by observed 'fireballs' in 
the galaxy RB199 in the Coma cluster by |Yoshida et al.| ( |2008 1. 
This galaxy shows a highly disturbed morphology that indicates 
a galaxy-galaxy merger. Interestingly, the galaxy also features 
filaments and bright blue knots at the ends of these filaments, 
which strongly emit Ha and ultraviolet radiation, indicating on- 
going star formation. Yoshid aet al.| ( |2"00 8) concluded, that ram- 
pressure stripping caused these knots to form. The masses of 
these knots are ~ 10 6 - 1O 7 M0, corresponding to the faintest 
dwarf spheroidals in the Virgo cluster (ibid.) and agreeing well 
with the masses found in the gas knots in our simulations, dis- 
cussed in Sect. 3.4 |Yoshida et al.| ( [200 8 ) found a typical SFR 
of around 10 3 Mq yr" 1 , which is lower than in our simulation. 
But because the knots of RB199 consist mainly of young stars 
formed there, having depleted the gas reservoir, they should be in 
a later evolutional state than in our simulations. Because RB199 
is in a post-starburst phase, introduced by the merger event, an- 
other possibility is that this galaxy already consumed a lot of 
its gas before the onset of ram-pressure stripping, which would 
explain the discrepancy in the observed and simulated star for- 
mation rates. 

Moreover, the filaments, consisting of older stars, are possibly 
left behind during the stripping, because they do not feel the ram 
pressure and are attracted by the gravitation of RB199. This can 
be well explained by our simulations, where we also found a 
fraction of stars falling back to the disc. 

We aim in the near future to compare observations based on our 
own data. In this project, we will use deep narrow- and broad- 
band imaging with VLT/FORS2 to trace in-situ star formation in 
the gas wakes of two disc galaxies in the Virgo cluster. 

5. Discussion and conclusion 



We carried out numerical studies on the effect of ram-pressure 
stripping and especially the influence of a galaxy's bulge. 
Thirty-one different high-resolution N-body/SPH simulations 
were conducted. The size of the galaxy bulge was varied 
and we used two different ICM densities of 10~ 28 gcrrr 3 and 
10~ 27 gcrrr 3 . In each case, the relative velocity of the galaxies 
was lOOOkms 1 and the temperature of the ICM was 10 7 K. 

- The size and mass of the bulge of the model galaxies were 
varied. While the total mass of all galaxies remained the same, 
mass was shifted from the dark-matter halo to the bulge in the 
initial conditions. Hence, a deeper potential well was formed. 
This induced a higher gas mass concentration in the centre of 
the galaxy, which increased the star formation. 

- After 2 Gyr of evolution in isolation, the spiral structure 
of the ISM developed into concentric rings in the simulation. 
When ram pressure was acting on a galaxy, however, the ISM 
was compressed, which by itself caused an enhancement of the 
star formation rate. Moreover, the ISM was stripped and the 
galaxy can lose most of its gas up to 90% if there was high ram 
pressure. In this case of high ram pressure acting, the process 
was basically completed within 200 Myr. Owing to the steeper 
potential well of the same galaxy containing a bulge, less gas 
was stripped as stated by the Gunn & Gott criterion (Gunn 

1972| ). As confirmed by previous investigations, the 



|& Gott 



inclination angle does not have a big influence on the amount 
of stripped gas. Only if the inclination was nearly edge-on the 
fraction of gas in the wake of the galaxy was much lower. A 
difference in the amount of stripped gas for distinct inclination 
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angles can be found when lower ram pressure is acting on the 
galaxies. 

- Star formation was enhanced when ram pressure was 
acting on a galaxy. In the considered ram-pressure scenarios, 
the star formation rate was up to four times higher than for the 
same galaxy evolving in vacuum. The higher the ram pressure, 
the more gas was stripped and therefore fewer stars were 
formed in the disc. On the other hand, many more dense star 
forming gas knots were present in the wake, and the higher 
the external pressure on these knots, the more stars were formed. 

- When the inclination angle of a galaxy was modified, the 
star formation rate did not change significantly when low ram 
pressure acted on a galaxy. Still, the highest SFR could be found 
when the galaxy was flying face-on through the ICM. In the 
high-density ICM, increasing the inclination angle resulted in 
a slight decrease of the star formation rate. Furthermore, in the 
edge-on case the star formation rate was lower at the beginning, 
but higher at the end of the simulation. 

- The mass fraction of newly formed stars in the wake 
compared to the total mass was much higher for stronger ram 
pressure. Again, because less gas was stripped if the bulge mass 
was higher, fewer stars were formed in the wake when a big 
bulge was present. In agreement with previous results, less mass 
of newly formed stars was present in the wake after 500 Myr 
of evolution if the bulge size of a galaxy was increased, and 
obviously also the mass in the wake was higher when the ram 
pressure was increased. 

- Different disc gas mass fractions do have a direct influence 
on the star formation rate, which decreases with a lower 
disc gas mass fraction. Interestingly, the enhancement due to 
ram-pressure was less pronounced for larger gas mass fractions. 
Apparently, because the gas is less dense when the gas mass 
fraction is low, ram-pressure can compress the gas to bring 
it into a star- forming regime. On the other hand, when the 
gas mass fraction and hence the density was already high, the 
compression effect was smaller. 

- Our investigations show that for high ram pressure up to 
25% of the newly formed stars were formed in the wake of 
the galaxy, whereas with low ram pressure, only a marginal 
amount of around 1% of all newly formed stars were present 
in the wake. When the ICM density was high, gas knots were 
formed from the stripped material in the wake. These high dense 
structures contain a cool core and hence many new stars were 
formed in these self-gravitating gas knots. With evolving time 
the gas knots are losing mass due to star formation processes. 
Although the gas knots were heated by stellar feedback, and 
gas flows out via stellar winds, the gas knots remained intact. 
Obviously, the cooling was very efficient. Additionally, self- 
gravity and the external pressure from the ICM kept the gas 
knots dense. Hence, the mass loss was almost completely due to 
star formation. The stars in these knots were not subject to ram 
pressure and were attracted by the galaxy itself. Because the 
stars are collisionless, they fall through the disc and back again, 
oscillating. Therefore, stars were present also in front of the disc. 

- The presence of these gas knots corresponds well to 
observations of the 'fireballs' in the wake of RB199 in the 



Coma cluster by Yoshida et al. ( 2008 ) and the dense blue gas 
knots in the wake of IC 3418 in the Virgo cluster ( |Martin et al.| 



2005). In both cases, the most plausible mechanism to form 
these gas knots is ram-pressure stripping, also confirmed by the 
simulations we performed here. 
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Fig. 7. Surface density of the ISM (green) and the newly formed stars (isolines) are shown for a ram-pressure stripping scenario of 
model galaxy lnb in an ICM with a density of 1CT 27 gcuT 3 . Note the colourbars for both distributions. Four different time steps are 
shown. New stars in the wake are formed in the dense gas knots. After 200 Myr many stars are present also in front of the disc. The 
stars formed in the wake are gravitationally attracted by the disc and due to the collisionless dynamics, these stars from the wake 
are falling through the disc. Furthermore, with the initial onset of ram pressure, the gas disc is pushed back from the stellar disc, as 
can be seen in the timestep after 50 Myr of evolution in the ICM. 
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